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Cracking the dark matter code at the LHC 
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We consider the decay of a generic resonance to two visible particles and any number of invisible 
particles. We show that the shape of the invariant mass distribution of the two visible particles is 
sensitive to both the mass spectrum of the new particles, as well as the decay topology. We provide 
the analytical formulas describing the invariant mass shapes for the nine simplest topologies (with 
up to two invisible particles in the final state). Any such distribution can be simply categorized by 
its endpoint, peak location and curvature, which are typically sufficient to discriminate among the 
competing topologies. In each case, we list the effective mass parameters which can be measured 
by experiment. In certain cases, the invariant mass shape is sufficient to completely determine the 
new particle mass spectrum, including the overall mass scale. 

PACS numbers: 14.80.Ly,12.60.Jv,13.85.-t 



The dark matter problem and the mystery of the feeble 
neutrinos greatly motivate the ongoing LHC searches for 
new physics in channels with missing energy. Alas, at 
hadron colliders like the LHC, deciphering events with 
invisible particles in the final state is notoriously difficult. 

The problem is schematically illustrated in Fig. [TJ 
which depicts the generic decay of some new heavy res- 
onance A into N v visible particles Vi and N x "invisible" 
particles Xi (neutrinos or dark matter candidates) which 
leave no trace in the detector. A priori we have no way of 
knowing the underlying physics behind Fig. [TJ and thus 
we are missing the answers to some very basic questions: 
1) How many invisible particles are in the final state? 2) 
What are their masses? 3) What is the exact topology 
(i.e. Feynman diagram) of Fig. [TJ are there any interme- 
diate resonances, and if so, what are their masses? 

Historically, the topic of mass measurements has at- 
tracted the most attention in the literature (for a review, 
see [J). Unfortunately, virtually all proposed methods 
suffer from two drawbacks. First, one must typically as- 
sume the correct decay topology for Fig. [TJ including the 
correct number N x of invisible particles. If this guess is 
incorrect, the method does not apply. This motivates 
us to address the issue of the correct decay topology 
and number of invisibles concurrently with (perhaps even 
prior to) the more traditional question of mass measure- 
ments [2j. Second, most methods for mass measurements 
utilize kinematic endpoints, where the available statistics 
can be rather poor (in the sense that the most popu- 
lated bins are rarely near the kinematic endpoint). Here 
we shall instead concentrate on the region near the peak 
rather than the endpoint of the kinematic distribution. 
(The exact shape in the vicinityof the endpoint does con- 




FIG. 1: The generic decay topology under consideration. 



TABLE I: The number of inequivalent event topologies as a 
function of 1 < JV„ < 4 and 1 < N x < 5. 
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tain information about N x but is difficult to mea- 

sure precisely in the presence of backgrounds and detec- 
tor effects.) Our main result will be the derivation of the 
analytical formulas necessary to analyze the full shape of 
the invariant mass distributions of the visible particles in 
Fig. [TJ including the location of the peak. We shall then 



demonstrate how those results can be used to determine: 
1) the number of missing particles; 2) their masses; and 
3) the associated event topology. 

Our setup is as follows. We consider the generic de- 
cay from Fig. [JJ without any prior assumptions about the 
decay topology or the number of invisibles. As seen in 
Table [TJ the number of inequivalent decay topologies pro- 
liferates very quickly as we increase the number of par- 
ticles in the final state. Let us begin with the simplest 
and most challenging case of N v = 2, postponing N v > 2 
to a future study Q . According to Table [TJ there are 
2 topologies with N x = 1, shown in Fig. [2ja,b), and 7 
topologies with N x = 2, shown in Fig. GJc-i). Our main 
goal is to analyze and contrast the V1V2 invariant mass 
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FIG. 2: The nine N v = 2 topologies with iV x < 2. 

distribution 1 in each of those nine cases. 

The differential distribution of the invariant mass to = 
m VlV2 will be described by an analytical formula 



dN 

— =f(m;M A ,M Bi ,M Xji 



(1) 



which is a function of the unknown masses. Given the 
general formula ([1} for f(m), we can easily obtain its 
kinematic endpoint 

E = max{m} (2) 

and the location P of the peak of the /(to) distribution 

f(m = P) = max {/(to)} . (3) 

Let us also define the dimensionless derivative ratios 
to™ d n f(m) 



Rn — — 



/(to) dm n 



(4) 



By definition, Ri = 0, as long as /(to) is continuously 
differentiable at to = P, while R2 parameterizes the cur- 
vature of /(to) at to = P. 

The parameters E, P and R n are in principle all exper- 
imentally measurable from the distribution (JTJ) . Tradi- 
tional studies [3] have always concentrated on measuring 



1 We note that the resonance A is in general allowed to be produced 
fully inclusively, with an arbitrary number of additional visible or 
invisible particles recoiling against A in the event. This precludes 
us from using the $!t measurement, since it will be corrupted 
by the invisible recoils, which leaves us with m vlV2 as the only 
viable observable to study. The related combinatorial problem 
of partitioning the visiblcs in the event was addressed in 2, 6]. 



just the endpoint E, failing to utilize all of the available 
information encoded in the distribution /(to). The end- 
point approach gives a single measurement ([2]), which is 
clearly insufficient to determine the full spectrum of res- 
onances involved in the decay chain of Fig. [1] Here we 
propose to invoke the full shape ([lj in the analysis. We 
envision that in practice this will be done by performing 
unbinned maximum-likelihood fits of (JT|) to the observed 
data. In order to illustrate the power of the method here, 
it is sufficient to consider just the additional individual 
measurements of P and R 2 . Since they are obtained from 
the most populated bins near the peak, we can expect 
that they will be rather well measured. More impor- 
tantly, the additional information about P and R2 might 
be sufficient to completely determine the mass spectrum 
(see eqs. (|I3II4[) below). But first we need to present our 
results for (Q]|4]) in each of the nine cases in Fig. [2] 

The topology of Fig. \^a). For a three body decay to 
massless visible particles, one has 



/(to; M a , M x ) ~ to A 1 / 2 (to 2 , M 2 a , M\ 



(5) 



wher 



A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2xz. (6) 
In this case 
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where 



6 1+ (l + 3a 



a = 2M A M X /{M 2 A + M 2 ) 



(7) 

(8) 

(9) 
(10) 



Contrary to popular belief, one can now solve for both 
masses M A and M x , given two of the three measurements 
([7]|9|). For example, using the peak location P and the 
endpoint E, we find 



M A = 



E IP 2- 3(P/Ef 



2 \EV 1- 2(P/E) 2 



M ¥ = - - 
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(11) 
(12) 



Eqs. (|11I12D offer a new method of determining both M A 
and M x , which is a simpler alternative to the Mt2 kink 
method of [|[ , since here we do not rely on the |?t mea- 
surement at all, and do not require to reconstruct the 
decay chain on the other side of the event. 

In fact, one does not even need an endpoint measure- 
ment, since the peak location P and the curvature R2 
are sufficient for this purpose: 
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FIG. 3: Distribution of Ma and M x found by a maximum- 
likelihood fit to eq. ([5]) in f 0,000 pseudo-experiments with 100 
signal events (left) or 1000 signal events (right). The input 
study point has M A = 550 GeV and M x = 50 GeV. 
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1 / 2 FIG. 4: The topology disambiguation diagram. The different 

(14) color-coded regions delineate the range of values for R2 and 
P/E spanned by each decay topology from Fig. [2] 



To the best of our knowledge (|13I14[) represent the first 
and only method in the literature which is capable of de- 
termining the complete mass spectrum in a short SUSY- 
like decay chain, without relying on any kinematic end- 
point measurements. 

In order to get a rough idea of the precision of these 
mass determinations, in Fig.[3]on the left (right) we show 
the results from 10,000 pseudo-experiments with 100 
(1000) signal events each. In each pseudo-experiment, 
the two masses M A and M x are extracted from a 
maximum-likelihood fit of the simulated data to the full 
distribution ([5]). Fig.[3]shows that, as expected, the mass 
difference is measured quite well, at the level of ~ 1% 
with just 100 events. At the same time, the mass sum 
(or equivalently, the absolute mass scale) is also being 
determined, albeit less precisely: at the level of ~ 30% 
(~ 10%) with 100 (1000) events. 

The topology of Fig. W{b). Here one obtains the cele- 
brated triangular shape 



Ma, M XI and M X2 can be simultaneously determined 
from a fit of cq. (|18|) to the data, just like in Fig. [3] [5]. 

The topology of Fig. W(d). The invariant mass distri- 
bution of the visible particles V\ and V2 is not affected by 
the emission of invisible particles upstream and so this 
case is equivalent to the topology of Fig. [^a). The cor- 
responding results can be obtained from ([S][5]) with the 
substitution A — > B, since now the role of the parent 
resonance is played by the intermediate particle B. One 
would then be able to determine independently Mb and 
M X2 , while Ma and M Xl would remain unknown. 

The topology of Fig. .\E(e). Similarly, this case is equiva- 
lent to Fig.|5[b), with the substitutions A — s- B\ , B -> B 2 
and x ~* X2- Once again, the emission of the invisible 
particle xi upstream is not observable. The only mea- 
surable parameter in this case will be the endpoint E. 

The topology of Fig. Hf/j. We find 
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Unfortunately, the masses enter the shape ([15]) only 
through the combination (|16|) . which is the single effec- 
tive mass parameter accessible experimentally. 

The topology of Fig. \&(c). The shape is more conve- 
niently given in integral form, which is easy to code up: 
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The explicit formulas for P and (fT5)) will be shown in Q . 
The important point is that in principle all three masses 
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In this case, out of the 4 input masses entering the topol- 
ogy of Fig. one can measure three independent de- 
grees of freedom, e.g. M A /M B , M Xl /M B and Mj,-M 2 2 . 
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The topology of Fig. \B(g)- The shape is described by 
A B ^Jx(s,M^M^) (24) 

and it is easy to see that the results are obtained from (|20I 
[23]) with the substitution Ma O — -M X2 . In particular, 
the three measurable parameters in this case can be taken 
as M Xl /M B , M X2 /M B and M\ - Mf . 

TTie topology of Fig. \B(h). This is the "sandwich" 

3e is given by 

< m < e' 7 ' E, 



topology studied in [9| ■ The shape is given by 



rjm, 



f(m) ~ ^ (25) 
[ m In (Elm) , e r ' E < m < E, 
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2M Bl M B2 
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and 



E = [e*(M%- MlJ(M 2 B2 - M*J/(M Bl M B2 )] 1/2 , (27) 



P = 



Ee" 1 , T)<1; 
Ee 1 , rj > 1; 



P 2 = 



00, »? < 1; 

1, r?>l. 



(28) 



The distribution (f2~5j) exhibits a cusp at the non- 
differentiable point m = e~ ri E. In this case, there are 
5 mass inputs: Ma, M Bl , Mb 2 , M Xl and M X2 , but only 
two independent measurable parameters: rj and E. 

The topology of Fig. W( i). This is the "antler" topol- 
ogy which was studied in [10| for the symmetric case of 
M Bl = Mr ; and M X1 = M X2 . Here we generalize the 
result in [10( to arbitrary masses and find that f(m) is 
given by the same expression (|25p . only this time 
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2M Bl M B2 

E = [e«(M Bi -Ml)(M B2 -M* 2 )/(M Bl M B2 )] 1/2 (30) 



and identical expressions (|28f for P and i?2- Just like 
the case of Fig. [U[h), out of the 5 mass inputs, r\ and E 
are the only two measurable mass parameters. Table |TT] 



summarizes the final tally of input particle masses and 
independent measurable parameters for each topology. 

Each topology from Fig. [5] also maps onto a restricted 
region in the (B,2,P/E) plane, as shown in Fig. 2] (for 
convenience, instead of R2 € (0, 00), in the figure we plot 
^ tan" 1 i?2 £ (0, 1)). For example, the cyan circle at (1, 1) 
marks the prediction for the two topologies of Fig.[2jb,e), 



TABLE II: The number of mass inputs N m for each topology 
in Fig. [2] and the number of independent measurable param- 
eters N p in the definition of f(m). 
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while the magenta dot at (0.5, 0.37) and the magenta ver- 
tical line correspond to the two topologies of Fig. [2]Jh,i). 
The blue (red, green, black) points refer to the topolo- 
giesof Fig.EKd) (Fig.^g), Fig.df), Fig. He)). Fig.@ 
demonstrates that with the three measurements E, P 
and i?2, one can already begin to constrain qualitatively 
the allowed event topologies. A more quantitative com- 
parison of the invariant mass shapes would have to em- 
ploy a quantitative measure, such as the Kullback-Lciblcr 
distance [ll|. Following [12j . we revisit the antler topol- 
ogy example of Fig.^i) studied in [Toj]. Table UTT1 lists the 
number of events which are needed to disfavor the wrong 
event topologies by a factor of 1000. The table shows 
that, apart from the twin scenario of Fig. &h) which has 
the same invariant mass shape, all other alternatives will 
be ruled out with sufficient statistics. 



TABLE III: The number of events needed to disfavor the 
wrong topology by a factor of 1000, if the data comes from 
the antler topology of Fig. [2]4) for the study point consid- 
ered in pxj: M A = 1500 GeV, M Bl = Ms, = 730 GeV and 
M xi = M X2 = 100 GeV. 



Data 


(a,d) 


(b,e) 


(c) 


(f) 


(g) 


(h,i) 


(0 


698 


37 


96 


275 


698 


00 



[1] A. Barr and C. Lester, J. Phys. G 37, 123001 (2010). 

[2] Y. Bai and H. -C. Cheng, JHEP 1106, 021 (2011). 

[3] E. Byckling and K. Kajantie, "Particle Kinematics", 

John Wiley & Sons, 1973. 
[4] K. Agashe, D. Kim, D. G. E. Walker and L. Zhu, Phys. 

Rev. D 84, 055020 (2011); G. F. Giudice, B. Gripaios 

and R. Mahbubani, Phys. Rev. D 85, 075019 (2012). 
[5] W. Cho, D. Kim, K. Matchev and M. Park, to appear. 
[6] M. Blanke, D. Curtin and M. Perelstein, Phys. Rev. D 

82, 035020 (2010); A. Rajaraman and F. Yu, Phys. Lett. 

B 700, 126 (2011); P. Baringer, K. Kong, M. McCaskey 

and D. Noonan, JHEP 1110, 101 (2011). 
[7] I. Hinchliffe, F. Paige, M. Shapiro, J. Soderqvist and 



W. Yao, Phys. Rev. D 55, 5520 (1997); K. Matchev, 
F. Moortgat, L. Pape and M. Park, JHEP 0908, 104 

(2009) . 

[8] W. S. Cho, K. Choi, Y. Kim and C. Park, Phys. Rev. 

Lett. 100, 171801 (2008). 
[9] K. Agashe, D. Kim, M. Toharia and D. Walker, Phys. 

Rev. D 82, 015007 (2010). 
[10] T. Han, I.-W. Kim and J. Song, Phys. Lett. B 693, 575 

(2010) . 

[11] S. Kullback and R. Leibler, Ann. Math. Statist. 22, 79 
(1951). 

[12] C. Athanasiou, C. G. Lester, J. M. Smillie and B. R. Web- 
ber, JHEP 0608, 055 (2006). 



